Fourier expansions for the potentials of lattices of charge 
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In this note we apply the Poisson sum rule [T] to obtain 
formal expressions for the Fourier coefficients of the potential 
of a lattice of generalized charge. Each generalized charge is 
assumed to contribute to the potential a term which depends 
only on the vector displacement from the charge's position. 
The coefficients are explicitly calculated for two types of in- 
dividual charge potentials: Coulomb and Yukawa. For one 
and two-dimensional Coulomb lattices, we find that the non- 
zero frequency components decrease exponentially both with 
the distance to the distribution and with frequency. The 
exponential convergence of these expansions indicates that 
a truncated Fourier series will often provide both a simple 
and accurate analytic approximation for these potentials. In 
particular, this result may often be applied to justify use of 
the continuous charge distribution approximation for obser- 
vation points further from the distribution than the spacing 
between the charges. This result also serves to explain nu- 
merical observations which have been presented previously 
[21 [31 |3]. The resulting expansions for lattices of Yukawa 
type charges are seen to be closely related to those of their 
Coulomb analogs. Yukawa potentials are physically realized 
in linearly screened Coulomb systems [5], for example, and 
so are of direct physical significance. We end the note with 
a brief study of finite and disordered lattices. 

The general form of the Poisson sum rule allows for a sum 
over a periodic lattice to be replaced by an integral over the 
volume of the lattice. For an arbitrary individual particle 
potential function /, the rule states 
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Here, n is the dimension of the lattice, V is the n-dimensional 
volume per unit cell of the direct lattice, the vectors R are 
the direct lattice vectors of the distribution, and the vectors 
G are the reciprocal lattice vectors of the distribution [6]. 
Shifting the origin of integration in Eq.lJTJ) to the projected 
position of the observation point r immediately gives formal 
integral representations for the Fourier coefficients of the po- 
tential. 

For a lattice of Coulomb point charges of charge q, the 
individual potential function is given by 



(2) 



The evaluation of the right side of Eq. fTJ) for the general two- 
dimensional Coulomb lattice is carried out in the appendix. 
In each dimension the integrals may be carried out readily 
and we find the following expansions for the potentials of one, 
two, and three-dimensional Coulomb lattices, respectively. 
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Here, <E>o is the potential one would obtain through the ap- 
propriate continuous charge distribution approximation, d is 
the perpendicular distance to the distribution, and the func- 
tion Kq is a modified Bessel function of the second kind. For 
large arguments, the function Kq may be expanded as [7] 
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It follows that the non-zero frequency components will indeed 
be exponentially damped with both the frequency and the 
distance from the distribution for one and two-dimensional 
Coulomb lattices. The three-dimensional expansion Eq.©, 
while not exponentially convergent, also converges at a faster 
rate than the original Coulomb sum. In all three cases we 
observe that potential oscillations perpendicular to high den- 
sity lattice planes dominate those oscillations perpendicular 
to low density lattice planes. We note that the one and 
three-dimensional Coulomb results Eq.([3]) and Eq.© have 
been discussed previously and have been applied to obtain 
expressions for a lattice's Madelung constant [H1H]- The gen- 
eral two-dimensional result Eq.Q appears to have escaped 
prior notice, however. 

The integrations required to obtain the potentials of lat- 
tices of Yukawa type charges, with individual particle po- 
tentials of the form gexp[— jr]/r, may also be carried out, 
though the evaluations are slightly more difficult. In the 
context of linearly screened Coulomb systems, the parame- 
ter 7 is proportional to the square root of the concentration 
of screening counterions. The Fourier expansions obtained 
through application of Eq.fl]) for the potentials of one, two, 
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three-dimensional Yukawa lattices are 
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It is interesting to observe that the non-zero frequency coeffi- 
cients above may always be obtained from the corresponding 
unscreened coefficients by replacing |G| by (|G| 2 + 7 2 ) 1 / 2 . It 
is also interesting that the screened systems' zero-frequency 
coefficients always have a mathematical form similar to their 
non-zero frequency coefficients. The Coulomb expansions 
may be obtained from these results by taking the limit 7^0. 
Some care needs to be taken with the zero-frequency terms, 
however. 

The potential of a finite lattice may often be approx- 
imated by that due to an infinite lattice of charge. This 
approximation will be reasonable provided the observation 
point is near enough to the lattice to feel negligible edge 
effects. We may approximate the distance, D, where edge 
effects start to become significant by setting the magnitude 
of the dominant oscillatory term of the infinite lattice's po- 
tential equal to the one over the distance to the edge of the 
distribution. If the characteristic spacing between charges is 
taken to be one and the projected distance to the edge is L, 
we find for one and two-dimensional Coulomb systems, 

D~l gi. (10) 

Numerically, it has been observed that as the observation 
point moves away from the charge distribution, there is often 
a relatively sharp transition region of distances over which 
the edge effects start to become significant. A plot demon- 
strating this transition for a linear array of Coulomb charges 
is shown in Fig. ([I]). A fit to the linear region of this data 
shows that the slope is approximately equal to —6.281. This 
is quite close to the value of — 2tt predicted by Eq.Q for 
an infinite line of point charges. For Coulomb systems, the 
transition distance D appears to be well described by the 
logarithmic fit 

D « 0.89 + 0.27 log L, (11) 

for an equally spaced, one-dimensional array, and by 

D « 0.34 + 0.22 log L, (12) 

for a square lattice, two-dimensional array. 

The potentials of lattices slightly perturbed away from 
their equilibrium positions may also often be approximated 
by the potentials of perfect, infinite lattices. Such a per- 
turbed lattice may be considered to be the sum of a perfect 
lattice plus a series of physical dipoles. The exponentially 
decaying terms due to the perfect lattice will dominate those 
terms due to the dipoles over some range which depends 
upon the size of the displacements. Outside of this range 
the dipoles will dominate and the error will go as a power of 
one over the distance to the distribution. We conclude that 
for observation points near large, clean physical systems, the 
results obtained here for infinite, perfect lattices may often 
apply. 




Figure 1: A plot of the logarithm of the difference between 
3?c,i(d,0, 0) and §c,i(<2, 0,0) versus d for a 2 x 6400 charge 
linear array placed on the y-axis and centered at the origin. 



A Two-dimensional expansion 

From Eq.([T]), we have 
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In the above, Jn is a Bessel function of the first kind. 
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